ASYMPTOTIC EXPANSION FOR THE NORMAL IMPLIED 
VOLATILITY IN LOCAL VOLATILITY MODELS 
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Abstract. We study the dynamics of the normal impHed volatility in a local volatility model, 
using a small-time expansion in powers of maturity T. At leading order in this expansion, the 
asymptotics of the normal implied volatility is similar, up to a different definition of the moneyness, 
to that of the log-normal volatility. This relation is preserved also to order 0(T) in the small-time 
expansion, and differences with the log-normal case appear first at 0(T2). The results are illustrated 
on a few examples of local volatility models with analytical local volatility, finding generally good 
agreement with exact or numerical solutions. We point out that the asymptotic expansion can fail if 
applied naively for models with nonanalytical local volatility, for example which have discontinuous 
derivatives. Using perturbation theory methods, we show that the ATM normal implied volatility 
for such a model contains a term ~ VT, with a coefficient which is proportional with the jump of 
the derivative. 
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1. Introduction. Local volatility models have been introduced some time ago 
[UIJ], see [3] for an overview, in order to model the price evolution of a financial asset 
in a way consistent with the known European option prices on that asset. 

According to the Gyongy theorem 4 , there is a unique one-dimensional stock 
process which can reproduce a given terminal distribution at each time horizon. This 
process is usually written in a way reminiscent of a log-normal process 

dS{t) = S{t)aD{S{t),t)dW{t) + S{t)n{t)dt (1.1) 

where cr£)(5(t),t) is the so-called local volatility. In the most general case, it depends 
both on the stock price S{t), and explicitly on time t. The process (|1.1|) reduces to a 
simple log-normal evolution for the case of a constant local volatility. 

The local volatility and drift can be determined empirically from the observed 
market prices of options [U [5] using the Dupire equation 

dTC{K,T) - ^al{K,T)dlC{K,T)+fi{T){C{K,T)-dKC{K,T)) (1.2) 

The drift /i(T) is found from the time dependence of the forward price F{T) as 

^iiT) ^~dT log F{T). (1.3) 

In principle this relation allows the determination of the local volatility in terms 
of option prices with all strikes and maturities. However, in practice this leads to 
numerical instabilities, and an equivalent relation is used, which replaces the option 
prices with the log- normal volatility [5] . 

Many particular cases of local volatility ajj {S{t), t) have been considered in the lit- 
erature, such as for example the CEV model 6 , and the quadratic volatility model [7]- 

The log- normal assumption implicit in (II. ip . and the quotation of the implied 
volatility as log-normal volatility, are justified for equities, for which a negative asset 
value S{t) < is unphysical. For interest rates this is not a restriction anymore, and 
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in fact normal volatilities are commonly used in quoting market volatilities as yield 
volatilities. There are a few advantages to using a normal volatility language for rates, 
the most important of which is the fact that normal volatilities are more stable under 
shifts of the rates. Under certain conditions, calibrating a model to normal volatilities 
improves the calibration stability in time. 

In this paper we consider the local volatility model for situations where a normal 
process is more natural, and replace the evolution equation (jl.ip with 

dS{t) = aDiS{t),t)dW{t) + ^{t)dt (1.4) 

In Section 2 we derive the Dupire equation for this process, and convert it to an ex- 
pression for the local volatility in terms of normal implied volatilities. In Section 3 
we use asymptotic methods to derive an expansion in powers of time for the normal 
implied volatility. The leading term in this expansion is shown to have a very similar 
form to that well-known in the context of log-normal rates (the so-called Berestecky- 
Busca-Florent (BBF) formula [S]). The same result holds for the term linear in time, 
and differences with the log-normal case appear first at quadratic order in time. The 
general method is illustrated in Section 4 on a few examples of local volatility models. 
In Section 5 we consider the case of non-analytic local volatility functions, and show 
on an explicit example that the asymptotic expansion can fail in such cases. Using 
perturbation theory methods we prove that for local volatility models where the lo- 
cal volatility has a discontinuous derivative, the normal implied volatility contains a 
nonanalytic term ~ Vt which is proportional to the jump of the derivative. This 
signals a generic failure of the asymptotic expansion for such models, and underlines 
the need for care in their application. 

2. Dupire equation for the normal volatility. Consider the process St driven 
by the local volatility model with local volatility aD{St,t) 

dSt^aDiSt,t)dWt+^ltdt (2.1) 

We work in the risk-neutral measure, where the drift ^ ft ^ Dt is the difference 
between the interest rate rj, assumed to be deterministic, and the dividend rate Dt, 
which is assumed to be paid continuously. Note the form of the drift term, which is 
different from the one usually adopted in local volatility models with a log-normal- 
type evolution (jl.ip . This change is required by later convenience, and it will allow 
us to express cr^) in terms of the implied normal volatility. 

This different choice affects also the relation between the drift /it and the forward 
stock price Ft, which reads 

/,(T) = dTF{T) . (2.2) 

Note the difference with the log- normal case Eq. (|1.3p . 

We will consider European options (calls and puts) on the asset St- The undis- 
counted price of a call with strike K and maturity T is given by the usual formula 

C{K,T)^¥.[{St-K)+\ (2.3) 

The option prices can be quoted either in terms of log-normal and normal implied 
volatilities. We will consider here the normal implied volatility aN{K,T), which is 
defined through the price of a call option assuming normally distributed terminal 
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stock price St — -Ft + o-nVTX, with X ^ N{0, 1) a normally distributed random 
variable. The result for the call price is [5] 

CMiK,T) = (Ft - K)n(^^^^) + -^e-'^^-^-'^'l^^'^-^^aMVT (2.4) 

^ unvT ' V27r 

We would like to find the forward equation satisfied by the prices of the call 
options C {K, T) under the process (|2.ip . This will give the generalization of the 
well-known Dupire equation to this case. We start by quoting the result 

drCiK, T) = \al{K, T)dlC{K, T) - fi{T)dKC{K, T) (2.5) 

which has a very similar form to the usual Dupire equation with the obvious substi- 
tution K^a'j-,{K,T) ajj{K,T). In addition, the form of the drift term is different, 
due to the different choice for the drift term in (|2.ip . 

Proof. The proof of (|2.5p proceeds in close analogy with the proof of the usual 
Dupire equation (|1.2p . see e.g. [3]- Start with the Fokker-Planck (FP) equation for 
the pdf of S{T) starting with the initial condition (p{S, 0) = S{S) 

dTv{S,T) = ]^dl[al{S)^{S,T)] ~ dsWMS,T)\ (2.6) 

The call price is expressed in terms of the pdf (^(5, T) as 

C{K,T)^ / dS{S - K)+^{S,T) (2.7) 
Jo 

Taking a derivative with respect to time and using the FP equation (|2.6p we find 

/•oo 

dTC{K,T)^ dx{x - K)dTip{x,T) (2.8) 
Jk 



K 

oo 



dx{x - K){-d,{f,{TMx, T)) + \dl{alip{x, T))} 



The first integral is equal to 



dx{x-K)d^{pL{T)Lp{x,T))^-n{T) dxip{x,T) (2.9) 
K Jk 



and the second one is 



/"OO 

/ dx{x-K){-dl{al^{x,T))=<7l{K,TMK,T) (2.10) 

JK 

where we integrated by parts a few times, and assumed that the pdf <p(x, T) vanishes 
sufficiently fast as x oo. 

Putting everything together we obtain for 



1 f°° 
drCiK, T) = -ol {K, TMK, T) + ^i{T) dx^{x, T) 

^ JK 



K 

\al {K, T)dlC{K, T) - ^i{T)dKC{K, T) 
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where in the second line we used dKC{K,T) = ~ dx(p{x,T). This completes the 
derivation of the Dupire equation (12.51) for the process (j2.ip . 

The relation (I2.5P allows the determination of the local volatility (Jd{S,T) from 
the market prices of the options with different strikes and maturities. As mentioned, 
this determination is not very stable numerically, due to the appearance of 0/0 insta- 
bilities for short times and near ATM strikes. For this reason an alternative approach 
is preferable, which replaces the option price with the respective implied volatililty. 

In the following we will derive such a relation in terms of the implied normal 
volatility. We follow closely the derivation in [3^ , adapting it to the problem at hand. 

We start by expressing the Bachelier option price in terms of the new independent 
variables 



y^K^Fr, w ^ a%T (2.11) 



We have 



Cn{v,w) = ~yN{-^) + e-^'/(2«,)^^ (2.12) 

Setting the call price C{K, T) equal to the Bachelier price imposes a special functional 
dependence of w on y and time 

C[K, T) = Civ + Ft,T)^Cn {y, w{y, T) = [aUy + Ft, T)T]) (2.13) 

The derivatives with respect to K, T in the Dupire formula can be replaced 
with derivatives with respect to y^w of the Bachelier price CNiy,w), allowing for 
y— dependence in the w variable due to the smile 

dTC{K,T) = {dTw)K{d^C{y,w))y + {dTy)K{dyC{y,w))^ (2.14) 

= (^^^ ^dioCNiy,w) - iiTdyCNiy,w) 

dKC{K,T) = dyC{y,w) = dyCN{y,w) + ^d.^CN{y,w) (2.15) 
d],C{K, T) = dlCN + 2^d^y^CN + ^d^CN + (^fdlCN (2.16) 

The higher derivatives of CN{y,w) can be simplified with the help of the relations 

d^CN = 2d^CN (2.17) 

dywCw — ——BwCn (2-18) 
" w 

^l^N = (|j2 - ^)^-^^ (2.19) 

Using these identities, the Dupire equation gives the following expression for the 
local volatility as a function of the normal volatility 

al[K ^y + Fr,T) = idw/dT)rc=y+^, + Mdw/dy) 

The numerator of this expression can be written even simpler as the total time deriva- 
tive of u;(2/, T): ^w{y + FT,T) = {dw/dT) 

K=y+FT + fJ^ridw/dy). This result extends 
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the well-known formula for the log-normal case (see e.g. Eq.(l.lO) in 3 ) to the case 
of the normal implied volatility. In the log-normal case an additional term is present 
in the denominator, equal to —jQidyw)"^ [3]. Converting from w to the normal implied 
volatility aN{y,T) we find the more explicit result 

^2 ^ T) + TdTaliv, T) + tirTdyaliy, T) 

^-'^^)dy^N{y) + ^{dyaN{y)? + \Tdlal{y)~^{dyal{ym 

_ a%iy, T) + Tdraliy, T) + tiTT{dyal{y, T)) 

(1 - ^)^y^^^y))^ + T<JM{y)dl<jM{y)] ^ ■ ^ 

It is instructive to compare these results with the corresponding results for the 
log-normal case. The corresponding expression reads (note that the natural variable 
for this case is the log-strike x = \og{K/FT), instead oi y = K — Ft) 

^2 T)^ a|g(x, T) + Tdralsi^, T) + ^iTT{^^als{^, T)) 

" ' {l^^^d^OBs{xW+TaBs{x)dlaBs{x)-\alsT^{d^aBs{xW 

(2.22) 

We observe that the 0{T) terms are the same, but the 0{T^) term present here 
is missing in the normal case Eq. p.2ip . Apart from this, the two expressions are 
identical, up to the different definitions of the independent variables x and y, as 
pointed out above. 

We will neglect the possibility of explicit time dependence of the local volatility, 
and consider only the case of a time-homogeneous local volatility, which depends on 
5*4 alone, such as in the CEV models and the shifted log-normal model to be discussed 
below. 

In the short time limit T — )■ the equation (I2.21[) simplifies by dropping the terms 
proportional to T, and setting Ft — >■ -fo on the left-hand side. The equation can be 
solved in closed form and the asymptotic normal implied volatility is given by 

I \ y K - Fo 

Jo aoiz) Jfo ao{L) 

This result is the analog of a corresponding asymptotic expression obtained in the log- 
normal case in [5] . In this limit the equation (|2.22p can be solved with the well-known 
result given by the BBF formula 

( ^ ^ \og{K/FT) 

<JBS.o{x) = .X rf, = dL ■ (2-2^) 

Jo <To{z) J Ft Lao(L) 

This expresses the fact that, in the asymptotic short-time limit, the log-normal implied 
volatility is the harmonic average of the local volatility. 

For strikes close to the ATM point \y\ ^ Fq, the integral in the denominator 
of (j2.23p can be approximated by the rectangular rule, which gives a very simple 
approximation for the normal implied volatility near the ATM point 

aN{K) ^ aD{\(K + Fq)) - y^{K - F^fal{]^{K + F^)) (—1^)" |^^^ (2.25) 

with K an undetermined point in the {Fq, K) interval. The error estimate corresponds 
to the Newton-Cotes formula of degree 2 [S]. 



6 



V. Costeanu and D. Pirjol 



This extends the famihar Hagan- Woodward approximation [101112) for the imphed 
volatihty in the log-normal case to the normal case, and allows the use of intuition 
familar from the former case also to the latter. For instance, the "local volatility skew 
is twice the implied volatility skew" statement holds for both cases, with the obvious 
correspondence between the local volatilities in the two cases. 

The asymptotic result for the normal implied volatility cro(y) and its derivatives 
are related to the local volatility <JD{y) at the point y = K — Fq = Q as. 

ao(0) = a,,(0), a^,(0) = iai,(0) (2.26) 



3 ' 6ai3(0) 
'(0) = ia-(0)- ^^^(0)^^(0) , (a^(0)) 



2az3(0) 4a2,(0) 

These relations can be easily obtained by taking the derivatives of the BBF formula 
(|2:23)) at the ATM point y 0. 

3. Asymptotics of the normal volatility. In this section we consider the 
inverse problem of determining the normal implied volatility cr^^y, T) for a given local 
volatility model, with a given local volatility a£){y,T). This problem is equivalent to 
the that of solving the evolution equation (|2.ip for the local volatility model. There 
are several approaches to this problem in the literature. 

• The method of asymptotic expansion in time T of the Dupire equation. This 
was already hinted at in the previous Section, where we derived the leading 
term of this expansion (|2.24p . This method can be extended to higher orders 
in T, see [HII^. 

• The heat kernel expansion [161 [13 HI]. This method uses the expansion 
of the Green's function for the Fokker-Planck equation (|2.6p . exploiting its 
similarity to the heat equation. This allows the use of asymptotic methods 
developed for the parabolic partial differential equation with space-dependent 
coefficients [201 [II]- 

The simplest local volatility model corresponds to a time-homogeneous local 
volatility cr£)(S'(t)), which depends only on S{t), but not explicitly on time. We 
will restrict ourselves to this case in the following. 

As mentioned above, the time-dependence of the forward Ft introduces time de- 
pendence even in this case. The drift term in the evolution equation can be eliminated 
by defining the forward asset price for a maturity T 

F{t,T) = S{t) + dsfi{s) (3.1) 

This follows the driftless process 

dF{t,T) ^ aDiF{t,T) ~ dsn{s),t)dWit) (3.2) 

This shows that, even if the local volatility does not depend explicitly on time, the 
evolution of F{t,T) has explicit time-dependence introduced through the drift term 
[22]. 

In the following we will use as starting point the equation (jl.4l) and construct 
the solution as an expansion in powers of time T. There are two ways to solve the 
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equation (|2.2ip for (jf^{y,T): working at fixed strike K, or at fixed y ^ K — Ft- 
In both cases the time dependence is made explicit by expanding in powers in T. 
Ultimately, we would like to find the expansion of the implied volatility as a function 
of strike 

aN{K, T) = a^{K) + ai{K)T + a2{K)T^ + ■■■ (3.3) 

Working at fixed y one would have to deal with additional time dependence on the 
left-hand side of (|2.21l) introduced through Ft in aoiK = y + Ft) which will have 
to be expanded too. In addition, the form of the final result expressed as function of 
strike would have the form a^iK, T) — ^ - ai{y + Ft)T^, and would mix orders in T. 

Both these inconveniences are avoided by working at fixed K. We will thus treat 
Eq. (|2.2ip as a differential equation in K, at each order in T. This requires that 
we make explicit the time dependence in the factor y ^ K ~ Ft appearing in the 
denominator of (|2.2ip . The simplest way to do this is to redefine y = K — Fq, and 

replace y ^ y— {Ft — Fq) = y — hqT — ^fJ-iT"^ • in the denominator of (I2.2ip . The 

resulting time dependence is also made explicit by expanding in T. Here we expanded 
also the drift jiT in a power series in T 

/iT = /io + /iir + M2'r2 + .-. (3.4) 

corresponding to the expansion of the forward 

Ft^Fo+ [ dtfit = Fo+fioT+l^iiT^ + --- (3.5) 

The resulting expansion of the expression on the right-hand side of ()2.21|) in 
powers of T has a remarkably simple form to all orders in T 



Recall that in this equation and below we use the definition y = K—Fq. For simplicity, 
we denoted dy derivatives with respect to y as primes. The function Hj{y) depends 
only on CFo{y), cri(y), • ■ ■ , o'j-iiy), but not on aj (y). It depends also on the drift terms 
with k < j. 

The functional form (j3.6p implies that the expansion terms crj(j/) can be deter- 
mined recursively, starting with cro(y) which is given by the BBF formula 

t \ y K -Fa 

^o(y) = p = rK (^-^^ 

Jo (Td(z) JFo aoiL) 

Assuming that all Uk{y) with k < j are known (and thus the function Hj{y) is 
known) , the coefficient aj (y) can be found by solving the differential equation obtained 
by equating the terms of 0{T^) on the both sides of (|3.6p 
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The solution of the hnear differential equation (|3.8p can be found by the method 
of the variation of constants 



y 

where Gj{y) is determined by the equation 

2?/ /cro(y)V'+^ 



a,{y) = a,{y)[^y^'cM (3.9) 



(^)-Ci(,)+i7,(,) = (3.10) 



The solution is 

CM - C0) - f 1(1 - {^y'H.M (3.") 

= C, (0) -if /' J, H, {z) . 

^ Jo 0-£l(z)(75(z) 

The integration constant Cj{0) is determined by the condition that <Tj{y) does not 
have a singularity at y = 0, or equivalently K = Fq. This requires that Cj{0) = 0. 

The final form of the solution is obtained by putting together the two factors, 
and is 



^i{y)=My)(^^Y^'\- r dz — HJz)] (3.12) 

The value of the j— th coefficient at y = depends only on the lower order coeffi- 
cients, and does not require an integration. To see this, let's examine the contributions 
of the different terms in the square bracket in the expansion p.6p at the point y = 0. 
The first term vanishes, and the other two cancel among each other. This gives a 
relation for aj (0) in terms of the lower order coefficients 

a,(0) = -ao(0)^^^ (3.13) 

Recall that due to the definition y = K — Fq adopted in this section, the y = point 
corresponds to the strike K — Fq, which coincides with the ATM point K = Ft only 
if the drift vanishes. 

This solves the recursion problem for <yj{y) in terms of the <Jk{y) with k < j. The 
only remaining problem is to find the function Hj{y). This can be done by expanding 
the expression (jl.4p . and requires only algebraic manipulations. 

3.1. The solution for (Ji{y). Here we illustrate the general method outlined 
above on the example of the leading 0{T) correction to the BBF formula. The 
inhomogeneous term Hi (y) is 

Substituting this into p. lip we get the integrals 

Jq <to{z) - za'{z) Jo ao(z)-zaQ(z) 
= \ogia{y) - ya'M + log(c7o(0)) = - log (—^M^) (3.16) 

Vcro(0)CTLi(y)/ 



+1 r i-y 
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and 



h-fd/:^{i-^)^[ (3.17) 

Combining everything gives for the first order (hnear in time) correction to the 
BBF formula 

\ 2 Vcrz5(y)cro(0)/ Jo ^cro(z) (Td[z)J J 
Expressed as a function of strike, the 0{T) correction to the normal implied volatility 



Y ^ 2 Vcrz5(t/)cro(0)/ 



IS 



a,{K) 



(3.19) 



Of course, this correction is well-known and has been derived in |17j in the context 
of the log-normal implied volatility, using a representation in terms of the process 
for the forward stock price. As shown, at order 0{T) the asymptotic expansion of 
the log-normal and normal implied volatilities are related by the simple replacement 
of the log-strike variable x — \og{K/ Sq) with the variable y — K — Sq. We have 
been unable to find an explicit result for the drift term in the literature, apart from 
Ref. [22], see Eq. (2.6) in this paper. Note however that the second term in p.l9p is 
different from Eq. (2.6) in [22] which has l/aQ{z) — l/aj){z) under the integral. 

The y = value of the first subleading correction p.lSp can be obtained from 
p.l3p . and depends only on the ATM local volatility and its derivatives 

'^i(O) = \<^U0yoi0) + Jmo^(^z,(0) - ao(0)) (3.20) 

= ^anm2aoi0)aU0)-M0)r]- 

The ATM normal implied volatility up to 0{T) is given by 

a^j^\K = Ft) = ao{K = Fo + y.^T) + a,{K = Fo + ^ioT)T (3.21) 
= ao{Fo) + (ai(Fo) + Aiofiol^o))?^ + 0{T^) . 

The absence of a drift contribution to ai {K — So,T) implies the following relation, 
true for any local volatility model: the price of a call option with strike K = Sq in 
the presence of a constant drift fi is equal to the price of the same call option in the 
absence of the drift, plus a known correction term ^/iT, up to terms quadratic in time 

C{K = So, T- /i) = C(i^ = 5*0, T; /X = 0) + i^^T + 0[T^) (3.22) 



The higher derivatives of the 0{T) correction at the K = point are given by 
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(up to terms proportional to the drift /ig) 

'^i(O) - ^{<Tl{Oyi;iO)+aniOy^iO)a'^{0) - ^^(O))^) (3.23) 
<^'m - 2^(l2al,(0)a-'(0) - AanmmO)? + 48(ai5(0))V^(0) - 13^^^ 



(3.24) 

3.2. The solution for 172(2/). The second order coefficient can be computed 
using (|3.12p and is given by 



(T2{y) = (tq{v) 



0-0(2/) 



dz 



2(7 D{z)al{z) 



H2{z) 



where the inhoniogeneous term in the equation lor the 0{T^) term is 



(75 (y (To 



-2Ny 



(3.25) 



(3.26) 



1 r 

7V4 



2Nyi^^y + aoiyK{yf 



My) 
CTo(y) 



2Nyi^)(^)' + a,iy)a'^{y)+aoiy)a'liy)} 
0-0(2/) ^0-0(2/)/ J 



+H^''\y) 



where N = 1 — jpfa) dy(To(jj) = t^tI^- The term H^^^y) contains the dependence on 
drift, and is equal to 



^liy) 



2/io 



K(2/)^i(2/))' 



' \ 2 



CTO 



o^oc^o 


22// 




_ 7V2 


N ' 





2^0^^ (2/) 2 1 /cr^(2/) 



iV3ao(2/) ''°iVnao(2/) 



V (Tn f y) / 



Mo 



(3.27) 
2y_(j'Q fai' 

iV2 



(-)' 



7VV(7o(2/) 



N\aoiy) 



The expression lor H2{y) can be simplified using the equation satisfied by cri(2/) 

|(^)' + 4^+i^i(y)=0 (3.28) 
N\cro(y)J aniy) 



We obtain 



^2(2/) =3 



0-0(2/) 



(<^o(2/)cro(2/))^ ^o(2/)^i(y) +Hit^)(^y-^ 



4iV4 



iV2 



(3.29) 



which gives the following result lor the second order correction to the normal implied 
volatility 



^2(2/) 



z dz\ 



a|,(z)K(z)]2 anizWiiz) 



1 



.2(7 D{z)a'^{z) 



8a4(z) 



2(7l{z) 



2aD{z)(7l{z) 



-H 
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Its value at K ~ Sq can be obtained from Eq. p.l3p and does not require an integra- 
tion. Neglecting the drift terms, it is given by 

<y2iK = So) = + ^ao(0)K'(Or + l^'liO) . (3.31) 

This can be expressed only in terms of the local volatility function (JD{y) using 
Eqs. (ESS), (ESI), (15:^ for cri(O) and <i(0). 

4. Examples. In this Section we apply the asymptotic method of Sec. [3] to a 
few particular cases, in order to test the numerical convergence of the solution. We 
compare the asymptotic expansion either to exact solutions, or to a numerical solution 
of the Dupire equation obtained by solving it numerically in Mathematica [11) . 

4.1. Shifted log-normal local volatility model. Consider the local volatility 
model with shifted log-normal dynamics, and a constant drift term 

dS{t) = (do + 2hS{t))dW{t) + fidt (4.1) 

This equation can be solved in closed form, and the solution reads 

S(t) = [So + ^)e^bww-2bh -^ + ^f dse^b(w,-w.)^2b\t-s) 3) 
zo Zb Jq 

The first term is log-normally distributed, while the second term, proportional to the 
drift, has a more complicated distribution. Keeping only the first term, the call price 
with strike K and maturity T can be expressed in terms of the familiar Black-Scholes 
price CBs{K,F,aBs,T) 

C{K, So,T)^CBs{K = K+^,F = So + ^,aBs = 2b, T) (4.3) 

The call price with nonzero drift has a more complicated expression, and is considered 
in Appendix IB] 

The leading asymptotic term for the normal implied volatility is obtained from 

K - Ft K - Ft 

ao{K) = -J, — = 2b (4.4) 

JFt a-a+2bL ° ao+2bSo 

The process (j4.1[) is invariant under the simultaneous shifts — > + (5, ctq — > (Tq ~ 2&5, 
so for simplicity we can choose Sq — 0. The results for nonzero Sq can be obtained 
by replacing K ^ K — Sq. 

Expanding in powers oi y = K — Sq one finds the leading order coefficient of the 
asymptotic expansion 



62 , 63 



ao{K)^ao + by-^y' + —^y^ + 0{y^) (4.5) 



where ao = ao + 2bSo. 

The correction of 0{T) is (Ti(y), which can be obtained from p.lSp . Its power 
expansion around the y = point is 



12 V. Costeanu and D. Pirjol 



O.Olp ' ' ' ^ 0.01 




0.01 0.02 0.03 0.04 0.05 0.01 0.02 0.03 0.04 0.05 

Fig. 4.1. The normal implied volatility af^{K) — 2bSo for the shifted log-normal model com- 
paring the exact solution (black) with the asymptotic solution (7q(K) (red dashed), and the solution 
including 0{T) corrections (solid blue). Model parameters: So = 3%,(to = 0.8%, = 0.1, and 
T = 10 (left), T = 30 (right). 



The ATM second order correction is obtained from p.l3p . Neglecting the drift, 
its expression is 



I , . , . 1 



1 



^2(0) = --ao(0)ff2(0) = ~-My) 3^ - 5a,{y)a'^{y) - MvWHy) + {cToiyK{y)f 

1 



20 



40 



aob^ ■ (4.7) 



Collecting ah terms, the total ATM normal implied volatility in the (driftless) 
shifted lognormal model is given by an expansion in b^T, whose first three terms are 

aN{K = So,T)^ao{l- \\?T + i^b^T^ + 0{{\?Tf)) . (4.8) 

This can be compared with the exact solution (|4.3p . which gives for the ATM 
normal volatility 

^AL_^ = g \N{h^) N{-h^)\ (4.9) 
Expanding in powers of time we get 



6 40^ ' 336^ ' 3456 

(4.10) 



T^™(r) = ao(l - -\?T + —{\?Tf {\?Tf + {\?Tf + 0{\?Tf ) 



The first three terms agree with the coefficients from the asymptotic expansion ob- 
tained above. 

To get a sense for the numerical accuracy of the expansion, we show in Figure 
14. li the leading 0(T°) term a^iy) term in the asymptotic expansion (red dashed line), 
together with the first subleading correction ctq (y) + To\ (y) (blue curve) , compared 
with the exact solution (black curve). The agreement of the asymptotic expansion 
with the exact solution is very good, already at 0{T). 

We show in Table 1 numerical results from the asymptotic expansion of the model, 
compared with the exact solution. The model parameters have been chosen as cto = 
3%, & = 0.2. The three columns show the error in the ATM normal implied volatility, 
with respect to the exact result. At zeroth order in the T expansion, the ATM implied 
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Table 4.1 

Results for the ATM normal implied volatility in the shifted log-normal model, compared with 
the exact result, as a function of maturity T. The parameters of the model are ctq = 3%, b = 0.2. 



T 










1 


0.04 


0.0199% 


-0.0001% 


0.0000% 


2 


0.08 


0.0395% 


-0.0005% 


0.0000% 


5 


0.2 


0.0971% 


-0.0029% 


0.0001% 


10 


0.4 


0.1885% 


-0.0115% 


0.0005% 


20 


0.8 


0.3562% 


-0.0438% 


0.0042% 


30 


1.2 


0.5058% 


-0.0942% 


0.0138% 



normal vol is (To = 3%, but at large times it deviates from this value. From Table ITTl 
this deviation can be seen to be as large as —0.5% for T — 30. 

Next we consider also the case with nonzero drift /i. This can be included in 
the expressions for the coefficients of the asymptotic expansion as discussed in the 
previous section, see Eq. (I3.19p . In particular, the ATM normal implied volatility, 
working to order O(T^), is given by (|3.21l) . which for the case considered here, reads 

<7n[K = So + fiT) =ao + i-h^ao + by)T + 0{T^) (4.11) 

6 

where the b^T term comes from the expansion of the leading order term tro (y) around 

y = o. 

We can check explicitly that this is reproduced by the exact solution of the model 
with constant drift, which is derived in the Appendix [BJ In Eq. (IB.35P we obtained 
the first two terms of the expansion of the call price with K ~ So in powers of the 
drift fj,. 

The model (|4.ip is related to the model in the Appendix |B] by a simple mapping 

t = 46V (4.12) 
Wt = 2bWr (4.13) 
^t = Sr + ^ (4.14) 

, = ^ (4-15) 
Then if xt satisfies the process dxt = xtdWt + jidt, then St will follow the process 

dSr = (cto + 2bSr)dWr + lJLdT . 

In Appendix |B] it is shown that the price of an option with strike equal to the 
initial asset value a;o is given by 

C(i^ = a:o,.To,t) = a:oErf(^) +iMt + 0(/i') (4.16) 
The corresponding result for the model (14. ip is 

C[K ^So,So,T) = ^!I±f^Erf (^) + i^T + 0(^^) . (4.17) 



From this result we can determine the normal implied volatility in the model 
(|4.1I) at the point K = Sq. This can be found by comparing (|4.17p with the Bachelier 
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formula with CTAr(ii: = 5o,r) = ao{l + ciT + 0{T'^)), which gives ci = -ifo^. This is in 
agreement with the result Eq. (I4.1ip for (Ti{K — So,T) obtained from the asymptotic 
expansion and confirms the absence of a term proportional to fi in (Ji{K = Sq,T). 

4.2. Stochastic volatility inspired local volatility model. As a second ex- 
ample, we derive the short-time asymptotics of the implied normal volatility in a local 
volatility model inspired by a stochastic volatility model. Consider the model 

dS{t) = ,y{t)crD{S{t))dWi{t) (4.18) 
diy(t) = -/v{t)dW2{t) . 

where the two stochastic drivers have correlation corr {Wi{t), W2{t)) — p. The initial 
condition is S'(O) = Ft, v{0) = 1. 

We would like to use the methods of Sec.|3]to find the short time asymptotics of 
the normal volatility in this model. This can be done by relating first the model (|4.18|) 
to the SABR model [T2|, and then using the well-known asymptotic local volatility 
of the latter to find the equivalent local volatility of the model (j4.18p . We start by 
making the change of variable 

ziSt) = r ^ (4.19) 

J Ft '^D[y) 

which transforms the first equation (|4.18p into 

dz ^ i^{t)dWi{t) + dnit (4.20) 

This is identical with the evolution of the log-price xt = \og{St/Ft) in the log- 
normal SABR model [12], so we can take over the asymptotic solution of the SABR 
model for the short-time asymptotic limit |12[ 113] . The following one-dimensional 
process has the same terminal distribution of x(t) as the two-dimensional stochastic 
volatility model at leading order in a short-time expansion [13l [14] 

aMz)^^^z'~2pjz + l. (4.21) 

Higher order terms in the asymptotic expansion of the SABR model have been ob- 
tained in [niniiiHiiin]. 

The process for Zt can be converted into a process for the original asset price S{t) 
by an application of the Ito lemma 

dSit) = aD{S{t))a,ii{z{S{t))dW{t) (4.22) 

where the function z{S) is given in (|4.19p . 

We can use now the short-time asymptotics (j2.23p applied to this one-dimensional 
model to derive the leading asymptotics for the normal smile in the stochastic volatility 
model (|i?T5]l 

fj^T^^ ^-^^ ^-^T ^--Pr 

a^K, T) = ^ = ^m^) ('-'^^ 

J Ft ao{S)<y,a(z(S)) Jo a„ff(x) ^' ^ " 



where 



D{x) = log (4.24) 
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and 



z{K) = [ 



Ft (^D{y) ' 



(4.25) 



The result (|4.23|) gives the short-time asymptotics for the normal implied volatility of 
the stochastic volatility model (|4.18|) . 

We consider in the following the case of the normal SABR model with constant 
local volatility <td{S) — ctq: a-nd study the normal implied volatility of the local 
volatility model (I4.22p . The argument of the local volatility depends on the variable 
z{S{t), S{Q)) — l/(To{S{t) — S'(O)), such that we obtain the local volatility model 



with y = S{t) — S{0). Although the original stochastic volatility model (I4.18|) is well- 
defined only for non-positive correlation p < [15) . we will use the local volatility 
(|4.22p for both positive and negative values of p, effectively ignoring its origin as a 
stochastic volatility model. 

One can use now the results of Section [2] to derive the normal implied volatility 
of this model as an expansion in T. The leading order result is given in Eq. (j4.23p . 
and the first subleading correction in Eq. (|3.18p . 

We show in Figure 14.21 the results for the normal implied volatility in the model 
(I4.26P at leading order in T, and including the 0{T) correction cri(y). These are 
compared with an exact numerical solution. We note that for moderate maturities 
the agreement obtained by keeping only the first subleading term is satisfactory. 

5. Nonanalytic local volatility. The asymptotic expansion for the normal 
implied volatility in powers of T presented in Sections 2 and 3 is applicable only for 
analytic local volatility functions aoiS). If aoiS) is not analytic, the asymptotic 
expansion fails, and new terms which are non-analytic in time can appear. 

To illustrate this phenomenon, consider a model with local volatility which has a 
discontinuous derivative 



According to the asymptotic expansion formula, the implied volatility cro(y) (and 
all higher order terms) depends only on the values of the local volatility in the interval 
{Sq,K), but not on its values outside this interval. For this case, it implies that the 
implied volatility (JN{y) for y > 0(y < 0) should depend only on bji{bL). We show in 
Fig. 3 the exact solution for the implied volatility aN{y), obtained using a numerical 
solution in Mathematica (solid lines), and the leading asymptotic results (To(if), from 
the BBF formula (dashed lines). The three curves correspond to = —0.1 (red), 
— Q (blue), and ~ 0.1 (black), keeping 0.1 fixed. We observe that the 

implied volatility for y > Q depends strongly on the value of , which determines the 
local volatility outside the region (5*0, if). Although for b^ — 0.1 the exact result is 
in good agreement with the asymptotic expansion, the agreement becomes far worse 
for 6l = 0,-0.1. 

In order to understand better the poor convergence of the asymptotic expansion 
in the case b^ "^bji, we consider in more detail the particular case of the model (|5.ip 
with —b^ ~bji = b, for which an exact solution is known |23| . The exact solution for 




(4.26) 




y<0 
y>0 



(5.1) 
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0. 0083 
0. 0086 ■ 
0.0084 ■ 
0.0082 

0.008 ■ 
0. 0078 ■ 
0. 0076 ■ 
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Fig. 4.2. The asymptotic expansion for the normal implied volatility of the SABR-like model 
{4-26^ compared with the exact numerical solution (black curves) for two maturities: T = 10 (left), 
and T = 30 (right). The red dashed curves show the asymptotic solution ao{K), and the blue curve 
shows the asymptotic solution including the 0(T) correction. Model parameters: So = 3%, ctq = 
0.8%, 7 = 0.2. The correlation is p = 0% (upper plots), p = 30% (middle plots) and p = —30% 
(lower plots). 



the ATM implied volatility can be written as a sum of two terms (see the Appendix for 
a detailed proof). The first term coincides with the solution of the shifted log- normal 
model considered above with h = h^, and the second term gives the correction due to 
the discontinuity at ?/ < 



CT(0,i) =CTi(0,<) + cr2(0,t) 



with 



(72(0, 



1 



-fJoe 2 



(5.2) 

(5.3) 
(5.4) 



Expanding the exact implied volatility in powers of time, the first term ai{0,t) 
reproduces the result of the shifted log-normal model (|4.10l) corresponding to the local 
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Fig. 5.1. The normal implied volatility af^{K) for the model \5. 1\) with parameters S = 
3%, (TO = 0.8%, 6i = O.l(black), O(blue), —O.l(red), bn = 0.1 and T = 10. The solid lines show exact 
numerical solutions of the Dupire equation, while the dashed lines correspond to the asymptotic 
solution ao(K). 



volatility in the y > extrapolated to the entire region ajj (y) = ctq + 2by 

ai(0, t) = ao(l - + ^bh' - ^bh' + ■■■) (5.5) 
However, there is a second correction a2{0,t), which has the time expansion 

The term (T2(0, t) is not reproduced by the asymptotic expansion discussed in Section 
3. Even more surprising, it contains a term of 0{^/t), which is not allowed by the 
asymptotic expansion of the Dupire equation in volatility form. This term largely 
accounts for the large deviation from the asymptotic result in Figure 3. The absence 
of the term a2{0,t) and its nonanalytic dependence on time signal a failure of the 
asymptotic expansion for the model (jS.ip . 

This failure appears to be generic for nonanalytic local volatility functions. We 
prove next the following result: in any local volatility model with a local volatility 
function which has a jump in the derivative at the forward point 

M'^UFt))' = dKaUK = FT + e)- Ok^tHK = Ft - e) (5.7) 

the time expansion of the ATM normal implied volatility contains a term ^ \/T which 
is nonanalytic in time, and which is proportional to the jump of the derivative 

CT7v(Ft,T) - cjd{Ft) + ^^(TD{FT)Vf/^{cjl{FT))' + 0{T) (5.8) 

This anomalous term reproduces the 0{\/T) term in Eq. (|5.6|) for the model (j5.1|) . 
for which the jump of the derivative is Acr^(FT) = 45. 

The proof of this result makes use of perturbation theory for linear operators [Ml 
[25] . We would like to solve the Dupire equation (|2.5I) which we repeat here for 
convenience. This gives the price of a call option C{K,T) at time T, and reads 

drCiK, T) - {K)dlC{K, T) (5.9) 
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with the initial condition C{K, 0) — max{K—SQ, 0). For simplicity, we assumed a zero 
drift term = 0, and assumed the local volatility to be time-homogeneous. This 
equation is identical with the heat equation with position dependent conductivity 
aniK), and the initial condition C{K, 0) = ma.x{K — 5o, 0). 
The solution can be written formally in operator form as 

C'iK, T) = exp {l{K)T^ C{K, 0) = exp (^LqT + (L - Lo)^) C{K, 0) (5.10) 

where we denoted the differential operators 

LoiK) = lal{So)d], (5.11) 

L{K) = lal{K)d'j, (5.12) 

The solution with L = Lq is well-known, and can be written as a convolution of the 
initial condition with the heat kernel 

Go{x, t; 2/, 0) = ^ exp(-^^— ^) (5.13) 
V zvrr 

The solution of the unperturbed equation dTCo{K,T) — Lo{K)C{K,T) is given by 
the convolution 

/oo 
dyGoiy, y, 0)max(y - So, 0) (5.14) 

where we denoted for simplicity crD('S'o) = ctq. The integral can be computed in closed 
form, and the result is the well-known Bachelier formula (|2.4p . 

We will treat the difference AL = L — Lq as a perturbation, and write the solution 
(|5.10p as an expansion in powers of AL 

CiK,T) ^ CoiK,T) + CiiK,T) + C2{K,T) + ■ ■ ■ (5.15) 

The term of zeroth order Co{K,T) is identical to the Bachelier result. The term of 
first order in the perturbation is given by a triple integral 

/oo /•oo 
dy du dTGo{K,alT-u,alT)AL{u)Go{u,alT-y,Q)inai^{v-So,0) 
-OO J — OO J 

(5.16) 

The action of the operator AL(ii) can be expressed as 

2 ( \ 2 

AL{u)Go{u,alT-y,0) - ^rGo(^, agr; y, 0) (5.17) 



(^Diu)-'7l ( 1 {u-y) X 2 



2t 2agT2 



Go(u, CTqt;?/, 0) 



Let us examine in some detail the dependence on time T of the first order per- 
turbation Ci{K,T) for the ATM case K = Sq. For simplicity, shift the origin of the 
y axis such that 5*0 = 0. It is also useful to introduce rescaled variables 
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Table 5.1 

Contributions to the ATM implied volatility from the n-th order of perturbation in AL. 



n 


ATM implied vol 


1 




2 




3 




4 





in terms of which the heat kernel is expressed as 

Go(0, agT; u, a^r) = ^^Go(0, 1; u, A) (5.19) 
We have in terms of these variables 

Ci(0,T)=aoVT r dydu f dXGo{0, l;u,X) '^oi^'^oVT) - gg ^q) 

J-oo Jo '^0 

x(-^ + ^^^2^)Go(S,A;y,0)(27)+ 
A similar rescaling gives for the zeroth order term 

/oo 
dyGo{u,X;y,Om+ (5.21) 
-oo 

We observe that the overall scaling factor ^/T is the same for Co and Ci, and the 
only new time dependence in Ci appears from the factor involving the local volatility 
a^{uaoVT)- Assuming analyticity, Taylor expansion of this factor produces both 
integer and half-integer powers of time 

aliuaoVT) -al = aoVTuiaHO))' + ^{aoVT)^u\aUO)y' + ■■■ (5.22) 

At this point we stop and examine the structure of the time-dependent terms at 
higher orders in perturbation theory. The n-th order in perturbation theory C„ {K, T) 
will be given by a 2n + 1-dimensional integral, which includes n double integrals over 
intermediate space-time variables, plus one final integral over the initial condition. 
The time dependence will consist of an overall ctoV^ factor as in Cq.i (if, T) plus terms 
arising from the presence of n factors of the form (j5.22p . The resulting contributions 
at each order n in perturbation theory can be easily obtained and are shown in Table 

o 

It is clear that only the first order in perturbation theory can contribute a term 
of 0{\/T) to the ATM implied volatility. However, upon computing the integrals over 
u, y we get a vanishing result 

dyyj duGo{0,l;u,X)ui-—+^-^^)Goiu,X;y,0)^ J dyfe'^y- (3 ~ f) = . 

(5.23) 

This proves the absence of a 0{\/T) term in the ATM implied volatility for the case 
of an analytical local volatility. This is in agreement with the results of Sections 2,3, 
which considered only this case. 
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Fig. 5.2. The normal implied volatility afj{K) for the model h5. il) comparing the exact solution 
(black) with the asymptotic solution (Jo{K) (red dashed), and the solution including 0{T) corrections 
(solid blue). Model parameters: S = 3%,(Jo = 0.8%, T = 10. 6^ n are as shown in the plots. 



However, the situation is very different if we assume a nonanalytical focal vofatif- 
ity. Taking an expansion for ajj{uao^/T) witti a finear term of tire form 

^2 fn-,^ ./^-\ 2 ( ucToVrbji if-u>0 , , 

1^ uaoV-i ol if m < 

the integration over u, y does not vanish anymore, and the result is proportional with 
the discontinuity of the derivative bn — 5i 

CiiO,T)^aoVT(^^aoVTibR~bL)) (5.25) 

Using the relation C{0,T) = ^aNiO,T)VT, we obtain from this result the 0{VT) 
term in the ATM implied volatility. This concludes the proof of the relation ()5.8p . 

The explicit result for the normal implied volatility Eq. ()5.6p for the model (j5.ip 
contains only one non-analytic term ~ VT, but not other terms of similar form, e.g. 
rp3/2^ rp5/2^ _ ^ would be interesting to investigate whether this is a general result, 
or if it holds only in the specific model (|5.ip . 

Although the relation (|5.8|) was proved for the ATM normal volatility, a similar 
result holds also for the ATM log-normal volatility. This is related to the former by 
the exact relation (|A.2|) 

.»(0,r) = f.;|Erf(f-M^) (5.26) 

In the short-maturity limit T — > this reduces to the well-known relation (7^(0) = 
^tCTbs(O)- The relation between <tn{0,T) and assi^iT) is an analytical function of 
time, which implies that a term proportional to ^/T in the expansion of the ATM 
normal implied volatility will introduce such a contribution also in the ATM log- 
normal implied volatility. 

Another manifestation of the failure of the asymptotic expansion for non-analytic 
local volatility functions is the prediction of a discontinuity in the implied volatility. 
This can happen for example in the model (|5.ip at order 0(T) in the T expansion, 
provided that 7^ The correction of 0{T) to the normal implied volatility is 
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Fig. 5.3. Discontinuity of the asymptotic solution at 0{T) for the model \5. with \b]^ \ 7^ . 
As explained in text, the subleading correction cri{y) is discontinuous at y = 0. The black curve 
shows the exact (numerical) solution, the red dashed curve shows the asymptotic solution a-o{K), 
and the blue curve shows the asymptotic solution including 0{T) corrections. Model parameters: 
S = 3%, cro = 0.8%, T = 10, 6i = 0.1, 6_R = 0.2. 

given by Eq. p.l9|) . The jump of the subleading correction can be obtained from the 
power series expansion of cri(y) in the shifted log-normal model, and is 

ai(0 + e)-a,{0-e) = -^{bj, - bl)ao (5.27) 

This discontinuity is shown in Figure 3 for a specific choice of the model parameters. 
Such a discontinuity is clearly an artifact of the naive application of the asymptotic 
expansion, as the numerical solution is continuous everywhere. 

6. Conclusions. In this paper we extended the short time asymptotic expansion 
method to local volatility models described in terms of normal volatility, as opposed 
to log-normal volatility. A description in terms of normal implied volatility appears 
naturally in the context of interest rates, which can become negative in regimes of 
small interest rates. The Dupire equation can be formulated as a nonlinear equation 
for the implied normal volatility, which can be solved by an asymptotic expansion 
in powers of time T. This equation is similar (although in a different independent 
variable - the difference strike y = K—Ft as opposed to the log-strike x = \og{K/FT)) 
to the usual BBF equation [5], from which it differs only at 0{T'^). 

We present explicit solutions for the coefficients of the 0{T), 0{T^) terms in the 
small-time expansion of the normal implied volatility. The drift term in the 0{T) 
coefficient can be expressed as a simple integral over the zeroth order coefficient and 
the local volatility. We point out that the drift contribution vanishes at the ATM 
point K — So in the 0{T) term. This absence is verified on the explicit example of 
the shifted log-normal model with constant drift, for which an exact solution can be 
obtained. 

We studied the convergence of the asymptotic expansion on two examples of 
analytical local volatility: shifted log-normal model with and without drift, and a 
model inspired from stochastic volatility models for which the local volatility is the 
square root of a quadratic polynomial of strike. We found generally good agreement, 
even when keeping only the first two terms 0{1), 0{T) in the small-time expansion. 

The asymptotic expansion of the implied volatility in integer powers of T can fail 
if the local volatility is a nonanalytic function. Using perturbation theory techniques. 
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we show that in local volatility models where the local volatility has a discontinuous 
derivative at the ATM point, the ATM normal implied volatility contains nonanalytic 
dependence on time, proportional to \fT . Furthermore, the coefficient of this term 
is simply determined, and is proportional to the jump of the derivative at the non- 
analyticity point. Although these results have been proven for the normal implied 
volatility, similar results must hold also for the log-normal implied volatility, which is 
related to the former by an analytical function of time. 

Another manifestation of the failure of the usiial asymptotic expansion for such 
models is the presence of a jump discontinuity in the 0(r) contribution to the implied 
volatility. These observations show that the usual asymptotic expansion can fail for 
non-analytic local volatility functions, and has to be used with care in such situations. 
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Appendix A. Relation between the normal and log-normal implied 
volatilities. 

We quote here an exact relation between normal and log-normal implied volatil- 
ities valid for arbitrary maturity T, but for small deviations from the ATM point 
y = K — Ft <^ Ft- This takes the form of an expansion in y/Fx, and the first few 
terms are 



WN 



- 



-e 



^ 2w 2V2^Erf A / \FtJ 

We denoted here WN{y,T) = a]^\/T and w{y,T) ~ (TbsVT the square roots of the 
normal and log-normal variances, respectively. 

The ATM implied volatilities satisfy the relation 



aN[0,T) ^ Frd ) (^-2 



2^/2 

which reduces in the short-maturity limit T — > to the well-known relation ctjv (0) = 
Fro-BsiO)- 

Taking the derivative with respect to strike at the ATM point gives an exact 
relation between the normal and log-normal skews 



a^v(0,T) = FT[a'Bs{0,T) + ^^ETi(^ 



aBs{0,T)VT 



2%/2 

which reduces in the short time limit to 



(A.3) 



^5v(^t) = \<yBs{FT) + Fra'BsiFT) (A.4) 

The asymptotic results of Sec. 2 give an exact relation between the normal and 
log-normal implied volatilities, valid in the short-time asymptotic limit. We limit 
ourselves to the case of time-homogeneous local volatility. The relations ()2.24p and 
(I2.23P express both of them in terms of the same local volatility anix). Eliminating 
the local volatility gives a direct relation between the two types of implied volatilities. 

The starting point is Eq. (|2.23p . where we substitute the local volatility auiy) 
with its expression in terms of the log-normal volatility 

aD[x) = — ^ -r^ (A.5) 

with X = log(i^/^T)- 

The integral in the denominator of (|2.23|) can be performed by integration by 
parts 

dL — \— = [ — - [ ^ ^ daBS = ^;r^ (A.6) 

and we find the very simple result 

aN{K) = aBsiK) ^~^^ (A.7) 
log TV 
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Expanding around this point one can find relations among the skews and curva- 
tures (convexities) of the normal and log-normal ATM skews, valid in the short time 
limit. We quote here the relation among convexities 

(^n{Ft) = —^aBsiFr) + (t'bs{Ft) + Fra'^siFr) (A.8) 

Appendix B. Log-normal process with constant drift. 

As an explicit example of local volatility model of the type (|1.4|) with non-zero 
drift, we consider the simplest such process, with log-normal volatility and constant 
drift 

dxt = XtdWt + ndt (B.l) 

We will compute the normal implied volatility in this model, and compare its small 
time expansion with the asymptotic expansion derived in Sec. |31 As noted in Sec. SI 
this explicit example confirms the correctness of the 0{T) correction to the normal 
implied volatility p.l9|) . 

The stochastic differential equation (jB.ip can be integrated in closed form with 
the result 

Xt = xoeMWt -lt)+n f dseMWt -Ws-^{t-s)). (B.2) 

From this expression it is clear that Xt is restricted to the range Xt > 0, provided that 
fi > 0. In the following we will restrict ourselves to the case ^ > 0. 
The mean and variance of Xt are 



E[xt] ^ xo + ^^t , (B.3) 
vniixt) = E[xt] - E[xtf = xl{e* - 1) + 2xonie* - 1 - t) + 2fi^{e* - 1 - t - —) . 



The process for zt = log(xt) is mean-reverting around Zmid ~ log(2/x). This reads 
dzf — dWt + (/ie~^' — ^)dt, which is similar to the Ornstein-Uhlenbeck process, but 
with an exponential mean-reverting term. From Eq. (jB.3[) we see that the variance 
of Xt increases with t and is unbounded, which is different from the OU proces, for 
which it approaches a finite limit as t — ^ oo. 

We would like to compute the price of a call option on xt with strike K conditional 
on xq ai t — 

C{K,t)^E[{xt-K)+\xo,t^O] (B.4) 
It satisfies the Dupire equation 

dtC{K, t) ^K^dj,C{K, t) - txdKC{K, t) (B.5) 

with the boundary conditions 

C(if,0) = (.To-/0+ (B.6) 
C(0, t) = xo + Ait , dKC{K,t)\K=f)^ -I, \im dKC{K,t)={). (B.7) 
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We will solve the Dupire equation (jB.Sp using Laplace transform methods. Define 
the Laplace transform of the call price with respect to time 

/>oo 

y{x,s) = / dte'^^Cix^t) (B.8) 
Jo 

This satisfies the Dupire equation in Laplace transformed form 

^x'^y"{x) - fiy'{x) - sy{x) = X - xo , < x < xq (B.9) 

= , X > Xq 

with boundary conditions 

y(0,s) = ^ + 4, y'{0,s) = --, \un y'{x,s)=0. (B.IO) 

S S S x—>oo 

The solution of (jB.9p is given by the sum of a particular solution of the inhomo- 
geneous equation plus the most general solution of the homogeneous equation 

( ^^{x~xo-^) + c,h{x)+c,h{x) x<xo ^^^^^ 

' \ ^2/2(2;) X>Xq ^ ' 

where fi^2{x) are the solutions of the homogeneous equation. They are 

=a;^iiFi(-Ai,2A2;-^) (B.12) 

X 

/2(x) =a:^^Fi(-A2,2Ai;-^) (B.13) 

X 

with iFi{a,b, z) the confluent hypergeometric function, and A1.2 are functions of s 



Ai = -(l + VT+Si) >0, A2 = -(1- VTTS^) <0 (B.14) 

In the X > Xq region we kept only the second solution /2 (a;) since it is the only one 
which satisfies the condition y'{x) — > as a; — > 00. For large values of the argument 
X the asymptotics of the solutions is 

hix)^x^\ f2{x)~x^' (B.15) 

such that /i(x) is increasing while f2{x) decreases approaching as x ^ 00. 
The wronskian of the two solutions of the homogeneous equation fi^2{x) is 

h{x)r^{x) - f[{x)f2{x) = (A2 - Ai)exp(-^) = -VTT8^exp(-^) . (B.16) 

X X 

The boundary conditions at a; = (jB.lOp require that the constants ci_2 satisfy the 
conditions 

ci/i(0)+C2/2(0) =0 (B.17) 

ci/{(0)+C2/^(0)=0. (B.18) 

This will have nonzero solutions for ci^2 only if the wronskian of the two functions 
vanishes at x = 0. Using the explicit expression in (IB.16P this is seen to be indeed 
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the case. Only one of the two equations above is independent; we will choose the first 
one. 

The limit of the functions fi^2{x) at x = is 

Another condition on the constants ci. 2,(^2 follows from requiring continuity of 
y{x) at X — xq 

cifiixo) + 02/2(2:0) - ^2/2(2:0) = -4 (B-20) 
Finally, another condition is obtained from the normalization of the integrated 

pdf 

dxy"ix) / dxicif'^ix)+C2f^ix)) + / dxd2f'^{x) - - (B.21) 



which gives 



ci/{(xo) + 02/^(2:0) - ^2/2(2^0) = \ (B.22) 



Together with (jB.17p and (jB.20[) this equation fixes the constants Ci,2,(i2- These 
constants are found by solving the equations 

ci/i(0) + 02/2(0) = (B.23) 
ci/i(2;o) + C2/2(.To) - rf2/2(.2:o) - "4 (B-24) 

Ci/{(2:o) + C2/^(xo) ~ ^2/2(2:0) = - (B.25) 

s 



The solution is 



/2(2:o) + (m/s)/2(2;o) _ /2(2;o) + {p/s)f^{xo) „_^2/i^ 



ifiMf2{xo)-h{xo)f^{xo})s sVYTTs """"^^xo^ ^ ■ ^ 

^ hmf2{xo) + {ti/s)f^{xo)] 

' f2mmXo)f2{xo)-fl{xo)mXo))s ^ • ^ 

, ,0 ,^TO ^ir(Ai)r(2A2) 1/2(2^0) + Ws)f^ixo)] ^ 2^1 
A2r(A2)r(2Ai) sVT+SS ^^xj 

, _ /i(0)/2(2:o) ^ /i(2:o)/2(0) + [/i(0)/^(2;o) - /{(.to)/2(0)](a^/-s) 

/2(0)(/i(xo)/^(2:o)-/((2;o)/2(2;o))s ■ ^ ■ ^ 

This completes the solution of the Laplace transformed Dupire equation (|B.9p . The 
call price C{K, t) is obtained by taking the inverse Laplace transform of the solution 
y{x, s) given in Eq. (IB. lip with the coefficients ci,2, (^2 given in (|B.26p . (IB.27p . (IB.28p . 

We consider next the explicit result for y{xo,s), which is the Laplace transform 
over time of the "ATM" call price (the quotation marks are a reminder that this is 
the ATM point only for the case of zero drift; in the presence of the drift the ATM 
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point corresponds to x — Xq + /i/s) 

y{xo,s) = 4 + exp(2AVxo) ^ ^/^(.o)][/i(.„) - ^/2(x„)] (B.29) 



M /, , cxp(2^/a:;o)\ , 

+ — =— )/l(xo)/2(xo)) 



+.o^^^f^t^t(-A„ 2A2; tF,(-A2, 2A,; + • • .(B.30) 

SVl + 8s Xq Xq 

where the ehipses are proportional to factors of (2/i)^-'-+^* and (2/1)-'-+^-'-+^*. 

This is easily expanded in powers of the drift fi. Keeping only terms linear in fi 
we find 

yi-o, s) ^ -o^=L= (1 + ^^^^^ + Oi,^)) + 4(1 + (B.31) 

= xo^=i= f 1 - + 0{^^')) + 4(i + 7n^) (B.32) 



sVl + 8sV a;o2s V 52^2 2y/l + 8s 

The first term corresponds to the usual Black-Scholes result for the ATM call. Its 
inverse Laplace transform can be found easily and is 

1 ^ Vt 



Erf 



sVl + 8s V2V2 

Its expansion in powers of t contains -x/t , 1 / -s/i , ■ • • ■ 

There is a partial cancellation of the 0(/x) terms, such that we get the "ATM" 
call price (in Laplace transformed form) 

1 , M 



y{xo, s) = xo^== + -4 + 0(A.^ M^) (B.34) 
svl + 8s 2s^ 

The inverse Laplace transform of this result is easily computed, and we find the 
explicit expansion of the "ATM" call price in powers of the drift 

C{K ^xo,t)^ xoErf(^) + + 0{^?) (B.35) 
2v 2 2 



V2^ V 24 640 21504 

We verified that the same result is obtained also using perturbation methods for 
linear operators applied to the Dupire equation (|B.5I) . considering the drift term as a 
perturbation. 

Appendix C. Proof of the relation (|5.2|) . 

We present here the details for the exact solution of the local volatility model 
with the local volatility 

„ / N / (To - 26y , y < . . 

= 1 .o + 2fey, y>0 ^^-1) 

Define the new variable 
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which follows the process 



dz{t) = dWit) - b{z)dt , b{z) = I 5 



-6, z<0 
2 > 



The original variable y is expressed in terms of z as 



f§(e2^--l), z>0 



(C.3) 



(C.4) 



The fundamental solution of the process (jC.3|) is known exactly, see Eq. (5.14) in 



p{z,t;x,0) 



1 



/27rt 
1 



/27rt 



, (X~Z — htY . , 

exp(-^ ^7 -) + be 

exp(2&x-^^^4-^) + ^e''' 



di;exp( ) 



2t 

{v~bt)\ 
dvexjp{ ) 



It is easy now to compute prices of call options using the relation 

/•OO 

C{So,K,t)^ / du{y{u)~ K)p{u,t;z{So),0) 

Jz(K) 



x>0,z>0 

, x>0,z<0 
(C.5) 

(C.6) 



In particular, the ATM implied normal volatility in Eq. (|5.2p was obtained from the 
ATM call price according to 



C{So, So, t) 



1 



27r 



<^ATM 



{So)Vt 



(C.7) 
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